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In the present paper the Renormahzation Group (RG) method is adopted as a tool for 
a constructive analysis of the properties of the Frobenius-Perron Operator. The renormal- 
ization group reduction of a generic symplectic map in the case, where the unperturbed 
rotation frequency of the map is far from structural resonances driven by the kick pertur- 
bation has been performed in detail. It is further shown that if the unperturbed rotation 
frequency is close to a resonance, the reduced RG map of the Frobenius-Perron operator 
(or phase-space density propagator) is equivalent to a discrete Fokker-Planck equation for 
the renormalized distribution function. The RG method has been also applied to study the 
stochastic properties of the standard Chirikov- Taylor map. 

§1. Introduction 

Recursive maps represent a useful and powerful tool to model and to facilitate 
the understanding of the physical processes taking place in complex nonlinear sys- 
tems. In particular, they are widely used to study the various transition scenarios 
from regular to chaotic behaviour in nonlinear dynamical systems,^ ''^''^ to simu- 
late physical systems exhibiting anomalous diffusion,^ or to analyze the underlying 
dynamics in time series with 1/f noise in their power spectrum.'^''^ Iterative maps 
provide a convenient and effective method to investigate single-particle dynamics in 
accelerators and storage rings .^''^ 

The extremely complicated behaviour of specific trajectories in chaotic systems 
strongly suggests a probabilistic approach to the dynamics. Instead of tracing an 
individual trajectory in phase space, one employs a statistical mechanics approach by 
means of a distribution function of an ensemble of trajectories. The Frobenius-Perron 
operator of a phase-space density (distribution) function, which sometimes is called 
the Transfer Operator of that function or a phase-space density propagator, provides 
a tool for studying the dynamics of the iteration of the distribution function itself. 
The iterative map yields complete information of how the value of an individual 
phase-space point jumps around during successive iterations, so that one gets a 
good sense of the point dynamics but no sense of how iteration acts on densities 
with support on sets in phase space. The latter gap is filled by the Frobenius-Perron 
operator, which provides a rule to determine how the evolution of densities over 
repeated iterations is accomplished. 

In the present paper we adopt the Renormahzation Group (RG) method for a 
constructive analysis of the properties of the Frobenius-Perron Operator. The basic 
idea of the method is to absorb secular or divergent terms of the naive perturba- 
tion solution into renormalized integration constants (amplitudes). As a result one 
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obtains an evolution law embedded in an evolution equation the renormalized am- 
plitudes must satisfy, which describes the large-scale dynamics of the system. The 
latter is usually called the RG equation. The stages in the renormalization group 
reduction of a particular physical system are quite general and well defined, which 
makes the RG method universal and independent on the concrete details of the 
underlying dynamics and physical processes involved .^''^'^'^'^ 

The paper is organized as follows. The next paragraph serves as a reminder 
of the basic definition, derivation and properties of the Frobenius-Perron operator 
and provides the starting point for the subsequent exposition. In Paragraph 3, 
we work out in detail the renormalization group reduction of a generic symplectic 
map in the case, where the unperturbed rotation frequency of the map is far from 
structural resonances driven by the kick perturbation. Paragraph 4 deals with the 
resonance structure of a symplectic map. It is shown that in the case, where the 
unperturbed rotation frequency is close to a resonance, the reduced RG map of the 
Frobenius-Perron operator is equivalent to a discrete Fokker-Planck equation for the 
renormalized distribution function. In Paragraph 5 we apply the RG method to 
study the stochastic properties of the standard Chirikov- Taylor map, and (re) derive 
the diffusion coefficient in the quasi-linear approximation. In Paragraph 6 a few 
examples concerning the application of the results thus obtained are worked out in 
detail. Finally, Paragraph 7 is dedicated to conclusions and outlook. 

§2. The Frobenius-Perron Operator for the Henon Map 

The Henon map is defined by the following expression!^ 



where 



' - sm LO cos uj ' 



is the rotation matrix for one period of the map. In terms of accelerator physics 
application this is equivalent to one revolution along the accelerator lattice. The 
frequency uj and the parameter S 

UJ = 27TU, S= , (2-3) 

are related to the unperturbed betatron tune ly and to the strength of the sextupole 
(cubic nonlinearity) perturbation Aq. Here / is the length of the sextupole, is its 
location on the azimuth of the machine and R is the mean radius. 
The Henon map can be written as 

Z„+i=<z„+i= ( % 2 ), (2-4) 
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where TZ^ denotes the transposed of the matrix (|2-2|) . The Frobenius- Perron opera- 
can be calculated explicitly. We have: 

/„+i(x,p) = U/„(x,p) = j dCdri6{X-C)6{P-v + Se)fn{C,v) 

= fn{X,P + SX^). (2-5) 
Introducing the formal small parameter e and the action-angle variables 

x = V2Jcosa, p=—V2Jsma, (2'6) 

with ^ 

J = -(x^+p^), a = — arctan ^— ^ , (2-7) 

we write the Probenius-Perron operator represented by equation H2-5|l in the form 

fn+i{a + uj,J) = fn{x,p + eSx^). (2-8) 

§3. Renormalization Group Treatment of the Frobenius-Perron 

Operator 

The generalization of the Frobenius-Perron operator ()2-8() for a generic symplec- 
tic map with rotation is straightforward. We have 

fn+iia + u;,J) = fn{x,p + edxVN), (3-1) 

where V]sr{x) is a potential and dx denotes partial differentiation with respect to x. 
In the case of the Henon map for instance, the potential has the form 

VNix) = (3-2) 

Equation ()3-l|) can be written as 

fn+iia + u;,J) = e^(^-^^)^-/„(a, J). (3-3) 

Since the potential does not depend on the momentum variable p, 

U = {dMdp - {dpVN)dx = {dxVN)dp, (3-4) 

where Ly is the Liouvillean operator associated with Vn- Therefore, equation 1)3 •3() 
becomes 

fn+i{a + uj,J)=e'^^fnia,J). (3-5) 

We assume that the potential Vn, written in action-angle variables can be split as 
follows 

VN{a,J) = Vo{J) + V{a,J). (3-6) 
Respectively, the Liouvillean operator can be written as 

Ly = Lo + L, (3-7) 
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where 

Lo = -u;v{J)da, L = {daV)dj - idjV)da, (3-8) 

and 

^v{J) = (3-9) 

First of all, we consider the case, where the rotation frequency uj is away from 
nonlinear resonances driven by the potential V. Following the standard procedure 
of the RG method,'^''^ we seek a solution to equation H3-5|) by naive perturbation 
expansion 

oo 

/„(a,J) = ^eV^'=)(a,J), (3-10) 
fc=0 

where the unknown functions fjt\a, J) have to be determined order by order. The 
zero-order equation 

/i2i(a + ^,J) = /r(a,J), (3-11) 

has the obvious solution 

J) = e-"'^^«F(a, J) = F{a - nuj, J). (3-12) 

To this end F{a, J) is an arbitrary function of its arguments, and will be the subject 
of the renormalization group reduction in the sequel. 
The first-order equation can be written as follows 

/i'i(a + uj,J)- fl^\a, J) = LvF(a - nuj, J). (3-13) 
Standard but cumbersome algebra yields the solution to equation H3-13() in the form 

fl^\a, J) = (nLo + C^)F{a - nu, J), (3-14) 

where 

= {daVM - {dMda. (3-15) 
Furthermore, the potential V^{a, J) is defined according to the expression 

/ in \ V ( /V*™" 

VJa.J) = V,{a--.j). Vi(a.J)=^ ^ 

Some of the details of the calculation can be found in Appendix A. 
The second order equation is 

fi^^a + uj,J)- fi^^Ha, J) = Ly/«(a, J) + ^F{a - nuJ, J). (3-17) 

Since we are interested in the secular solution of equation (|3-17j) . we retain on its 
right-hand-side only terms that would yield a secular contribution. Thus, the second 
order equation giving rise to secular solution can be written as 

1^ 



/(5,(a + u;,J)-/f)(a,J) 



n + - ) Lq + nLLo + J)da 



F{a — nuj, J), 

(3-18) 
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where 



Q{uj,J) = ^ mcot ( — \dj{VmdjVm). 



in=l 



(3-19) 



Omitting the details of the calculation (presented in Appendix A), we can write the 
second-order solution as 



n 



Lg + n£^Lo + nQ{uj, J)da 



F{a — nuJ, J) + non secular terms. 



(3-20) 

To remove secular terms (proportional to n and n^) in the first-order (|3-14j) and 
the second-order solution H3-2U() . we define a renormalization group transformation 
F{a, J) — > F{a, J; n) by collecting all terms proportional to F{a — nu;, J) 



F[a — nw, J; n) 



1 + enLn + 



n 



Lg + nQda 



F{a — nu), J). 



(3-21) 



Solving perturbatively equation H3-21() for F{a — nuj, J) in terms of F{a — nuj, J; n), 
we obtain ^ 

F{a-nuJ,J) = (l - enLo + ■ ■ ■) F{a - nto, J;n). (3-22) 

Following ReferenceP'^ we define a discrete version of the RG equation by consid- 
ering the difference 

F{a — nu!, J; n + 1) — F{a — nu), J; n) 



= jeLg + e^ 

Substituting the expression for F{a — nto, J) in terms of F{a — nuj, J; n) from equa- 
tion (|3-22|) . we can eliminate secular terms up to 0(e^)- The result is 

F{a — nuj, J; n + 1) — F{a — nuJ, J; n) = 

(3-24) 

Equation ()3-24() is the RG equation. It describes the evolution of the distribution 
function on slow time scales in addition to the fast oscillations with a fundamental 
frequency lj. 

An important remark is in order at this point. Note that once the renormaliza- 
tion transformation has been performed, the second term in the second-order solution 
H3-2U() is eliminated as well. Combining it with the second (non secular) term in the 
first-order solution (|3-14() . we obtain 

eC^F + e^nCSoF = tZ^ (\ - enLg) F(n) + e^nZjL^F(n) = eZ^F{n). (3-25) 

To first order in the perturbation parameter e the renormalized solution to equa- 
tion ()3-5() can be written as 

/„(a,J) = (\^eZ^F{a-nuj,J-n), (3-26) 



n + 



F(a — nuj, J). 



(3-23) 



eLg + 




F{a — nu), J; n) 
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where the renormahzed "amphtude" F{a — noj, J; n) satisfies the RG equation (|3-24() . 
In the continuous hmit equation (|3-24p acquires the form 



dF{a — nto, J; n) 
dn 



F{a-nLO,J;n). (3-27) 



Provided Lq 7^ (in the case, where the potential Vn is not antisymmetric) the latter 
is a Fokker-Planck equation with the Fokker-Planck operator acting on the angle 
variable only. A relevant example of a cubic map will be considered in Paragraph 6. 

§4. Resonance Structure of a Symplectic Map 

The solution (|3-14|) to the first-order perturbation equation (|3-13|) was obtained 
under the assumption that the unperturbed betatron tune i' is sufficiently far from 
any structural nonlinear resonance of the form rnQi^ = 1, where mo is an integer. In 
the present paragraph, we assume that 

2tt 

oj = (jjQ + e6i + 6^62 + . . . , 000 = — • (4-1) 

mo 

Moreover, for the sake of simplicity, we assume that there are no higher angle- 
dependent harmonics in the Fourier spectrum of V{a, J) that would drive higher- 
order resonances of the form pniQV = p, where p is an integer. However, results can 
be generalized easily to take into account this case as well. 

Proceeding as in Paragraph 3, we write the zero-order solution as 

/(0)(a, J) = e-"-°^«F(a, J) = F{a - ncoo, J). (4-2) 

The first-order equation in the resonant case can be written as follows 

/i'^^i(a + a;o,J)-/«(a,J) = (Li+L)F(a-nu;o,J), (4-3) 

where 

Li = -6ida + Lo- (4-4) 

The solution to equation H4-3() is readily found in a standard manner, analogous to 
that already used in Paragraph 3 (see also Appendix A). The result is 

/«(a,J)= [n(Li + L^) +C]^(a-nu;o,J), (4-5) 

where 

L/J = {daVR)dj - {djVR)da, (4-6) 
is the resonant Liouvillean operator and 

VR{a,J)= ^ Vm{J)e'"''' = 2Vmo{J)cosmoa, (4-7) 
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is the resonant potential. Furthermore, 

C'^ = {daV!:)dj-{djV:,)da, 

where now the potential V^(a, J) is defined according to the expression 



VL{a,J) = Vi{a-'^,j)., 



Vi{a,J)= 



2i sin (ma;o/2) 



The second-order equation in the resonant case can be written as 



fi^l,{a + ujo, J) - fii'>{a, J) = -6idaf^'i,{a + cvo, J) 



P(2) 



(1) 



(4-8) 



(4-9) 



-^{^191 + 2525a)/i2i(a + LOO, J) + UA'\a, J) + ^F{a - ncoo, J). (4-10) 

Similar to Paragraph 3, we again retain terms on the right-hand-side of equation 
(|4-lUp that would yield secular contributions to the second-order solution. Thus the 
second-order solution can be written as 



F{a — nujo, J) 



H ^Li + Ijiij F{a — nuJo, J) + n'^'uj (j-'i + ^lij -^(o — nujQ, J) + non secular terms, 

^ (4-11) 

iR,da\ =LRda-daiR, (4-12) 



where 



is the commutator of the operators Jjr and da, and 

. .^ ^ ______ fmu;o' 



'{loo, J) = Yl "^^ot {-^)dj{VmdjV^). 



(4-13) 



m=l 



is the new nonlinear tune shift. Here the prime in the above sum implies that the 
term with m = mo is excluded from the sum. 

Repeating the steps that brought us along from equation (|3-21|) to equation 
(|3-24j) in the previous paragraph, we obtain the RG equation in the resonant case 

F{an,J;n + 1) - F{an,J; n) 



e Li +Li? +e 



[Q' -52)da + %\LR,da 



F{an,J;n), 



(4-14) 
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where an = a — nujQ. Analogously to Paragraph 3, to first order in the perturba- 
tion parameter e the renormalized solution of equation (|3-5|) in the resonant case is 
represented by the expression 

fn{a,J) = (l + tCl]F{a-nuJo,J-.n), (4-15) 



As in the non resonant case, the renormalization group transformation ensures 
the elimination of the last term on the right-hand-side of equation (|4-lip [combined 
with the last term on the right-hand-side of equation H4-5() ] up to O(e^). In the 
continuous limit the RG equation (|4-14|) is a Fokker-Planck equation with a Fokker- 
Planck operator acting on both action and angle variables. 

§5. Stochastic Properties of the Standard Map 

The present paragraph is dedicated to the analysis of the stochastic properties 
of the standard map. The Frobenius-Perron operator can be calculated explicitly 
and the Frobenius-Perron equation can be written aP 

A+i(a, J) = e-^^-^e^^^-^^Ua, J), (5-1) 

or 

e'^^fn+i{a, J) = e"*'^'^^Vn(a, J). (5-2) 
We introduce a new variable ^ according to the equation 

i='^^, J = Jo + Vet (5-3) 

where Jo is a fixed value of the action variable. The first case we will consider is the 
one, where Jo ^ 0. The Frobenius-Perron equation (|5-2|) can be rewritten as 

e(^«+^«)^"/„,+i(a,0 = e^^''--^^Ua,0- (5-4) 
The zero-order equation 

e^«^»/i°i(a,0 = /Si(a + Jo,e) = fi'\a,0, (5-5) 
has the obvious solution 

(a, = e-^'°^-F{a, = F{a - nJo, 0- (5-6) 
The first-order perturbation equation can be written as follows 

e'''>^^fi%{a,0-fi'Hci,0 = {-Cda + smad^)F{a-nJo,0. (5-7) 
The latter can be solved in a straightforward manner and the result is 



cos (a - Jo/2) „ 
2sm(Jo/2) 



F{a-nJo,^). (5-8) 
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Proceeding further, we write the second-order perturbation equation 



pJoda /■(2) _ f{2) _ _ Joda 

^ Jn+l Jn — ^ 



(5-9) 

Retaining terms on the right-hand-side of equation that would give rise to 

secular contribution, we cast the latter in an explicit form according to the expression 

e'°^''£l, - = e[n + ^) dlF{a - nJo, - nsinad^^daFia - nJo, 

_coUJo/2)^.^2a4F(a-nJo,0. (5-10) 
The secular solution of equation ()5-10() is 

(a, 6 = ^edlFia - n Jo, + n^^^^^-^d^CdaF{a - n Jo, 0- (5-11) 
2 2sm(Jo/zj 

Following out a renormalization procedure similar to the one performed in Para- 
graphs 3 and 4, we obtain the RG equation 

F{a - nJo,^;n + 1) - F{a - nJo,(.;n) = (^-Ve(,da + ^dl^^F{a - nJo,^;n). 

(5-12) 

Note that in the case of Jo = tt, the third term on the right-hand-side of equation 
(|5-1U[) vanishes and we are again left with the above expressions for the secular 
second-order solution ()5-ll|) and for the renormalization group map ()5-12|) (taken for 
Jo = vr). 

To complete the present paragraph, we consider the case of Jo = 0. Equation 
H5-3() implies a simple (non canonical) scaling of the canonical variables, assuming 
that J is a slow variable. Then, equation (|5-lj) can be rewritten as 

fn+i{a,0 = e-^<^^e^^^--^^Ua,i). (5-13) 

Note that now the perturbation parameter is -y/e rather than e. Using the Campbell- 
Baker-Hausdorff identity 

exp (^oA) exp (^/5b) = exp ^qA /3B ^ A, B +..^, (5-14) 

for any operators A and B and any parameters a and /3, we cast equation (|5-13|) in 
the form 

/„4.i(a,e) = exp (V^L5 + |C + ...)/„ (a, e). (5-15) 

Here 

Jjs = —£,da + sina5^, C = [sinad^,^da] = sinada — £, cosad^. (5T6) 
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Next we proceed with solving of equation (|5-15j) order by order. 
The zero-order solution is trivial to find: 



(5-17) 



where as before, F{a, ^) is a generic function of its arguments, playing the role of an 
integration constant. The first-order equation 



tfiK6-/4'^(«,e)=LsF(a,0, 



has the obvious solution 



/«(a,0=nLsF(a,0. 



The second-order equation can be written explicitly as 



Its solution is readily found to be 



1 



fP{a,0 = ^(n'il + nC)F{a,0. 



(5-18) 
(5-19) 

(5-20) 
(5-21) 



Performing the already familiar renormalization transformation in analogy to what 
has been done in the previous paragraphs, we obtain the RG equation 



Fia,^;n + l)-F{a,C;n)= ^eLs + -(U + c) F{a,^;n 



(5-22) 



Averaging the RG equation H5-22[) over the angle variable a, we obtain a discrete 
version of the diffusion equation 



FoiC; n + 1) - FoiC; n) = -dlFoi^; n 



(5-23) 



for the angle-independent part of the distribution function Fo(Ci^) iii quasi-linear 
approximation. Higher order contributions to the diffusion coefficient e^/2 can be 
obtained by direct analysis of the Frobenius-Perron operator.^ 

§6. Examples 



Let us now consider a few examples. In the non resonant case of the Henon map 
Lq = ) equation ()3-27() can be written in the form 



dF{a — nuj, J; n) 



where 



dn 



Qh{<^, J) 



i^nia, J)daF{a — nuj, J; n), 



3 cot 



+ COt^-jj 



(6-1) 
(6-2) 
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Equation (|6-1|) written in alternative form 

^^^g/'""^ = [-a; + QH{a, J)]daF{a, J; n), (6-3) 
describes regular motion with a frequency uj — Qh-, and effective Hamiltonian 



52 

ife//(J) = WJ- — 



Scot ( - I + cot j 



J2. (6-4) 



The next example we would like to consider is the case of a cubic map with a 
potential of the form 

Vm{x) = ^. (6-5) 
Taking into account equation HS-6|) . we also have 

Vq{J) = ^—, v{a,J) = —-cos2a + ^cos4a. (6-6) 

8 2 8 

Furthermore, the nonlinear tune shift can be expressed according to 

3CJ SC'^J^ 
ujy{J) = —^ f2c{oJ, J) = —^{8 cot UJ + cot 2uj). (6-7) 

Thus, equation (|3-'27() can be written as 

^^^Ipl = -5a,F(a, J; n) + ^dlF{a, J; n), (6-8) 

where 

3(^2 j2 

a;(w, J) = w + (8 cot a; + cot 2a;). (6-9) 

32 

Equation ()6-8() can be readily solved, yielding the result 

F{a, J;n) = Y, Fk{J] ^y^^'^-'^) e-^^^'l^'l^ . (6-10) 
k 

The latter indicates that the renormalized distribution function -F(a, J; n) rapidly 
relaxes towards the invariant density Fq{J). 

To complete the present paragraph, we study the resonance case for the Henon 
map. We assume that the unperturbed betatron tune is close to a third order 
resonance Svq = 1. For the operator Li and for the nonlinear tune shift f2'{uJo,J), 
we obtain 

Li = -5ida, f2'{u;o, J) = (6-11) 

o 

respectively. In the continuous limit equation ()4-14|) can be written as 

1^ = _ _ f2')dj + LjiF + ^ (tj + 2LiLr + L|) F, (6-12) 
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where 

L_R = — ^^^^[2Jsm {3a)dj + cos {3a)da], (6-13) 

the frequency co is represented by expression (|4-1|) up to O(e^), and the renormahzed 
distribution function F = F{a, J; n) is a function of the phase-space variables and 
the "time" n. 

Equation (|6-12|) is a Fokker-Planck equation describing the slow evolution of the 
phase-space density in the case where the rotation frequency of the Henon map is 
close to a third order resonance. 

§7. Concluding Remarks 

We have applied the Renormalization Group (RG) method to study the stochas- 
tic properties of the Probenius-Perron operator for a variety of symplectic maps. Af- 
ter a brief introduction and derivation of the Probenius-Perron operator for a generic 
symplectic map with rotation, the case, where the unperturbed rotation frequency of 
the map is far from structural resonances driven by the kick perturbation has been 
analyzed in detail. It has been shown that up to second order in the strength of the 
perturbation kick, the renormalized propagator for maps with nonlinear stabilization 

^Lq 7^ 0^ describes random wandering of the angle variable. Purther, the resonance 
structure of a symplectic map has been investigated. It has been shown that in 
the case, where the unperturbed rotation frequency is close to a resonance, the re- 
duced RG map of the Probenius-Perron operator (or reduced phase-space density 
propagator) is equivalent to a discrete Pokker-Planck equation for the renormalized 
distribution function. 

The RG method has been also applied to study the stochastic properties of the 
standard Chirikov- Taylor map. A nontrivial discrete analogue of the Pokker-Planck 
equation with a Pokker-Planck operator acting on both canonical variables has been 
obtained. The latter reduces to the well-known diffusion equation in quasi-linear 
approximation for the angle-independent part of the distribution function. 

It is worthwhile to mention that the procedure developed in the present paper 
(see Paragraph 4) can be applied with a slight generalization to study modulational 
effects in symplectic maps. 

Appendix 

Writing equation (|3-13|) in the form 

/(^i(a + u;,J)-/«(a,J) = (Lo + L)F(a-n^,J), (A-l) 

we notice that the right-hand-side will give rise to two kinds of terms. Let us first 
consider the equation 

ipn+i{a + uj,J)- (fn{a, J) = ^QF{a - nuj, J). (A-2) 
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Taking into account the commutativity between Lq and da, we can rewrite equation 
as 

e^^^^n+i{a, J) - (/^„(a, J) = e-^^^^UF{a, J). (A-3) 
It is straightforward to verify that the solution of the last equation is 

ipn{a, J) = ne-'^^^-toFia, J) = nh^Fia - nw, J). (A-4) 

Since the potential y(a, J), the arbitrary function F{a,J) and the first-order 
distribution function fn\a, J) are periodic in the angle variable a, we can represent 
them clS cl Fourier series in a 

V{a, J) = Y. Kn( J)e*"^ F{a, J) = Fs( J)e^^^ (A-5) 

/«(a,J) = j;4")(J)e^^^ (A.6) 

k 

We substitute the above expansions into both sides of the remainder equation 

ipn+i{a + uj,J)- V'n(a, J) = LF(a - nuj, J), (A-7) 
and after equating similar harmonics, we obtain 

G'^'^-'^e^'^- - = ^ [imVmFl,_,, - i{k - m)!^^^.^] e"*^^-'-)"-, (A-8) 

m 

where the primes indicate differentiation with respect to the action variable J. It is 
straightforward to verify that the solution of equation HA-8|) has the form 

= E 

m 

(A-9) 

Substituting back expression ()A-9|) into the expansion (|A-6|) for the function ^pn and 
rearranging terms, we obtain 

^„(a, J) = 

m,s 

Since equation l\A-l\\ is linear, its general solution can be represented as a sum of 
expressions ()A-4() and (|A-10() . To complete the derivation of the first-order solution 
H3-14() , we note that expression (|A-10|) represents the Fourier expansion of the second 
term in H3-14|) . provided the potential Kj(a, J) is given by expression (|3-16|) . 

Here, we briefly sketch the derivation of equation ()3-2U|l . The first and the last 
terms on the right-hand-side of equation ()3-18() can be treated in a way analogous 
to the treatment of equation HA-2() . Consider the solution of the equation 

'Fn+i{a + 10, J) - 'Fn{a, J) = nLLoF{a - nuj, J). (ATI) 



im 



T/ „—imuj/2 
'me 

2i sin {muj/2) 



FL 



i{k 



m] 



T^/ —imuj/2 

2i sin {muj/2) 



Fk 



-i{k—m)nu) 



T/ -imu}/2 

2i sin {muj /2) 



F' 



is^^^ —Fs 

2i sin [mu>/2) 



(ATO) 
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Using the representation ()A-5|) and (|A-6|) . we can write the solution of equation 
(|JaT|) as 

^„(a,J) = ^g(")(J)e^'=^ (A.12) 

k 

(n) 

It can be verified by direct substitution that the functions Q^, are given by the 
expression 

= E (^^'^"^ + ^fcm)e-*('^-™)"", (A-13) 

m 

where 

"TA p~imuj/2 jr/ p—imuj/2 

Mm = ^m-^ - - ^W^fc-m, (A-14) 

Bum = im . ,^. Wl_^ - i{k - rn)—^f^—Wk-rn. (A-15) 
Asm (mu)/2) Asm [mu)/ 2) 

Here 

W{a,J) = LoW{a,J). (A-16) 

Thus, the desired expression for the second term on the right-hand-side of equation 
(|3-2Ujl is readily obtained by taking into account the Fourier expansion of the solution 
of equation (|A-11|) . 
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